Universality of the Scaling Law for Ferroic Domains 
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We show how the periodicity of 180° domains as a function of crystal thickness scales with the 
thickness of the domain walls both for ferroelectric and for ferromagnetic materials. We derive an 
analytical expression for the universal scaling factor and use this to calculate domain wall thickness 
and gradient coefficients (exchange constants) in some ferroic materials. We then use these to discuss 
some of the wider implications for the physics of ferroelectric nano-devices and periodically poled 
photonic crystals. 
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The generic term "ferroic" designates crystalline ma- 
terials that are ordered either ferroelectrically, ferro- 
magnetically or ferroelastically (including also antifer- 
roic configurations). Ferroic materials usually display 
domains, that is, regions that are either ordered along 
different polar directions or along the same direction but 
with opposite polarity (180° domains). Kittel showed 
some 60 years ago that the width of 180° magnetic do- 
mains (w) is correlated to the thickness of a crystal in 
a very well defined manner: the square of the domain 
width (w) is directly proportional to the thickness of the 
crystal (d) Kittel's law was latter extended by Mitsui 
and Furuichi (1953) for ferroelectric materials [2j, and by 
Roytburd (1972) for epitaxially clamped ferroelastic ones 

a. 

Recently, Schilling et al. have shown that the con- 
stant of proportionality between w 2 and d is a defin- 
ing characteristic of the type of ferroic transition be- 
ing considered, with ferromagnets having generally big- 
ger domains than ferroelectrics for crystals of the same 
thickness. These experimental results were rationalised 
by Scott who observed that the differences between 
ferroelectric and ferromagnetic domain periodicity essen- 
tially disappeared once the domain wall thickness was in- 
corporated as a scaling factor. Mathematically this was 

2 

expressed as — M, where T is the thickness of the 
domain wall and d is the crystal thickness; since domain 
walls tend to be narrow (few unit cells) for for all ferro- 
electrics, and broader for ferromagnets (tens of nanome- 
ters), the dimensionless factor M ends up being pretty 
much the same for both. This is nicely illustrated in Fig- 
ure 1: the different characteristics for w 2 as a function 
of crystal thickness of ferroelectrics and ferromagnets fall 
into the same parent curve once the square of the domain 
width is scaled by the domain wall thickness T. 

In this Letter we discuss the nature of the dimension- 
less constant M, and derive a simple analytical expression 
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for it. We then use the result to calculate the thickness 
of domain walls in several ferroic materials, and discuss 
some of the wider implications for practical applications. 

We begin by writing the Landau thermodynamic po- 
tential describing a second-order phase transition in a 
ferroic. For simplicity's sake, we shall focus on the case 
of 180° walls in a uniaxial ferroic, so a single order pa- 
rameter will suffice. We call this order parameter Q, and 
it can be either the polarization or the magnetization, de- 
pending on the type of ferroic being considered. Across 
a domain wall there is a change of the sign of the order 
parameter. That means that there is necessarily a gradi- 
ent, whose associated energy must be incorporated into 
the thermodynamic potential. Assuming that Q points 
along the z direction and that the domain wall is per- 
pendicular to the x axis, the thermodynamic potential is 

aG =i« 2 +^+KS) 2 w 

For a mono-domain state, and also for the center of the 
domains, we can neglect the gradient term, and mini- 
mization then leads to the familiar result for the order 
parameter in the ferroic state: 

Ql = ~ (2) 

The second derivative of the free energy with respect to 
the order parameter is the stiffness. In the ferroic state, 
the result is x _1 = —2a or, relating this to the order 
parameter, 

xc 1 = ml (3) 

Here the term "stiffness" has different meanings de- 
pending on the ferroic context, being inverse permittivity 
for ferroelectrics and inverse susceptibility for magnets. 

The energy density of the domain wall is calculated by 
minimizing the energy difference between a mono-domain 
state and a state with one domain wall. That is, one has 
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FIG. 1: (Color online)Above:Square of the 180° domain width 
as a function of crystal thickness for some ferroics (data ex- 
tracted from refs. HEim). Below: When the square of 
the domain width is divided by the wall thickness, all data- 
points fall essentially on the same parent curve. The domain 
wall size used for the scaling has been extracted from refs. 
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to minimize the following 0,0: 



dx 
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Where Qq is given by eq. (2). Minimization of a with 
the boundary condition Q(x=±oo) = Qq leads to the 
solution of the profile of the order parameter across the 
domain wall Q(x)=Qo tanh(x/<5), where the characteris- 
tic thickness is 



(5) 




and the energy density of the domain wall is 



parameter to just one dimension we have implicitly dis- 
carded the possibility of Bloch walls and Neel walls. This 
is in principle wrong for magnetic materials (though not 
for ferroelectrics). However, the analysis of such walls 
in magnetic materials actually arrives at the same solu- 
tions [jj. Second, in limiting our Landau expansion to 
the Q 4 order, we are limiting ourselves to second order 
transitions, which is not the case for several important 
ferroelectrics such as PbTiC>3 and BaTiC>3, although in 
thin films of these two materials epitaxial clamping (re- 
gardless of strain) changes the transition from first to 
second order || Q ■ For the exact solution of the first 
order domain wall the reader is adviced to look at ref. 
[lfH . And third, we have neglected the elastic coupling 
to the lattice distortions (strain terms), which is partic- 
ularly important in ferroelectrics as they generally are 
also ferroelastic. However, the effect of strain can be in- 
corporated by a renormalization of the coefficients in the 
Landau expansion |fj, |8|, l2l 111! , so 0UT treatment is still 
be valid once the renormalized coefficients are used. 

Regarding the physical interpretation of eq[51 k repre- 
sents an "exchange" constant, as its energy contribution 
is proportional to the mismatch of spins/dipoles with re- 
spect to their neighbors (gradient term) whereas a and 
b represent the "anisotropy" contributions, as they indi- 
cate the strength of the alignment of the order parameter 
with respect to the crystallographic axes. Quite natu- 
rally, it follows that if the anisotropy terms are big, the 
domain walls will tend to be narrow so as to minimize 
the number of misaligned spins/dipoles, whereas if the 
exchange k is big the domain walls will tend to be wide 
so that the gradient is as small as possible. In magnets, 
the exchange interaction wins, whereas in ferroelectrics 
the opposite is true, hence the difference in domain wall 
thickness generally observed between the two types of 
ferroic 0,0- 

The exchange constant k is well characterized for 
most magnetic materials, but that is not the case for 
ferroelectrics|l2). a problem which has so far complicated 
analysis based on eqs. |5]and[|)] It is therefore interesting 
to write the energy density as a function of domain wall 
thickness, which removes the dependence on k: 

a=^bS = ~ 

This expression will be used later. 

We turn now to the relation between domain period- 
icity and thickness in a crystal slab cut perpendicular to 
the polar direction. The uncompensated dipoles/spins at 
the surface generate a large electrostatic/magnetostatic 
energy, which is reduced by creating domains of opposite 
polarity. The depolarization/demagnetization energy of 
the two crystal surfaces as a function of domain width is: 
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Several important simplifications have been made in 
the above treatment. First, in constraining the order 
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Where £(3) is Riemanns zeta function £(3) ~ 1.202 13J . 
Although the physical forces involved (electrostatic, mag- 
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netostatic) are different, the Maxwell equations for the 
energy are analogous, and thus the resulting expression 
for the surface energy ends up being much the same 
IH 0- S fl5L ITt^ : the difference between the two fer- 
roics is thus not contained in the shape of the equation, 
but only in the magnitudes involved: the order param- 
eter Q (polarization/magnetization) and the permittiv- 
ity/susceptibility x- 

The reduction in surface energy achieved by introduc- 
ing domains is of course partly offset by the energy cost of 
the domain walls, which is proportional to a, to the wall 
size (itself proportional to the crystal thickness d) and 
to the number density of domain walls (inversely propor- 
tional to the domain width w). Hence, F wa u — crd/w. 
Adding the two energy components and minimizing with 
respect to the domain width w leads to the standard re- 
sult 
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If we now substitute the order parameter Qq and the 
energy density a by their respective expressions from eqs. 
121 and [7| the final expression for the dimensionless factor 



is: 



M 




2.455, 



(10) 



The experimental observation that M is generally a 
number in the range 1-10 for any ferroic is thus explained: 
the result is always a numerical constant (~2.455) mod- 
ified by the square root of the susceptibility anisotropy. 
The appeal of this expression is not just in its simplic- 
ity and generality, but also in that it can actually be 
exploited for practical purposes. 

It has been a long standing and challenging problem 
to establish the thickness of domain walls in ferroics, and 
very specially in ferroelectrics, as the later tend to be 
very thin and hard to measure experimentally p"slll9l.l2Cj| . 
On the other hand, theoretical approaches based on phe- 
nomenological models suffered from the fact that the co- 
efficient k of the gradient term is hard to characterize 
experimentally. Eg uat ion 1101 dispenses the need to know 
such a coefficient. Measuring domain width and crystal 
thickness and knowing the dielectric constants of a mate- 
rial should be enough to estimate the domain wall width. 
So does it work? 

In figure 1 we have shown the square of the domain 
width as a function of crystal thickness for 180° domains 
in ferroelectric PbTiO^ and Rochelle Salt, and ferromag- 
netic Co. All of them can be analyzed with the present 
treatment, although in the case of the PTO films a correc- 
tion due to the effect of the substrate must be taken into 
account 0,EJ. The measured slope, the permittivities 
and the calculated thickness of the 180° domain walls are 
shown in table 1, next to previous values extracted from 
the literature. 
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TABLE I: Experimental slope of w 2 vs d for 180° domains 
in two ferroelectrics (slope=<5M), calculated thickness of the 
domain walls and calculated value of the exchange constant 
k. The Landau coefficients used in the calculations have been 
extracted from refs. and The domain wall thickness is 
compared with previous published estimates 



Our predicted value for 8 of the Rochelle Salt is 13A 
(and thus the wall thickness is T = 2<5=26A), compatible 
with the results of Mitsui and Furuichi (T=24-47A)0, 
and Zhirnov (<5=12-220A) 6j. As for the predicted value 
for the domain walls of PTO, once the effect of the STO 
substrate has been taken into account |l5l Il7| we obtain 
that 5=2. 45A or T=4.9A, in excellent agreement with the 
first principle calculations of Meyer and Vanderbilt[2l|. 
The combination of our model with experimental data 
thus agrees with previous estimates, and support the 
view that ferroelectric 180° domain walls are atomically 
sharp [HIH. 

The above equations apply to magnetic materials too. 
Fitting the data of Co to our model, the calculated do- 
main wall thickness is 20nm, which is somewhat thicker 
than previous theoretical estimates that yield a value of 
5-10nm, but thinner than the experimentally determined 
values of 46nm |24|. All in all, the results suggest that 
the method is quite robust for the analysis of different 
ferroics. 

Moreover, once the domain wall thickness has been 
determined, one can go back to eq. |S] and determine 
the value of the exchange constant (k) for the material, 
which is also a long standing challenge in ferroelectrics 
|12| . It is worth emphasizing here that this constant is 
of great importance in determining the performance of 
ferroelectric thin films and nanostructures ja, |26| , where 
the gradient term associated with the surface depolar- 
ization has a strong effect on the functional properties 
EHH H3. Our calculated values of k for PbTi0 3 

and Rochelle salt are included in table 1. 

Parenthetically, we note also that regular domains in 
ferroelectric crystals have important applications in pho- 
tonics, where they are used for frequency conversion 
through quasi phase matching . Presently, the reg- 
ularly spaced stripe domains are achieved through peri- 
odic poling, which has limitations due to the large coer- 
cive field -and some times finite conductivity- of some of 
the most important photonic crystals, such as LiNbOj, 
and KTiOPOi (KTP). Importantly also, artificially fab- 
ricated domain structures are generally not in thermo- 
dynamic equilibrium, and switchback can occur |34j . In 
theory it should be possible to achieve self-patterned and 
stable regular domains in ferroelectric photonic crystals 
by cutting them at the right thickness and preventing 
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charge screening upon cooling through the phase transi- 
tion, although in practice periodic poling is always likely 
to be required (e.g. to maximize registration). As a prac- 
tical example of reverse engineering, the measured value 
of the domain wall thickness for KTP is 3 A 19], while the 
longitudinal and transverse relative permittivities are 15 
and 11 respectively [35l |: accordingly, a periodic domain 
structure with a domain width of e.g. 5 microns would 
be most stable for a crystal ~3cm thick or, conversely, 
a 0.5mm crystal of KTP can have domains as small as 
0.7/im. This suggests that the known difficulty in sta- 
bilizing narrow domains in thick crystals is not due to 
intrinsic factors; indeed, experiments that make use of 
the depolarizing field have achieved self-patterned sub- 
micron domains in LiNbC>3 |3t| . 

In sum, the universal scaling law for ferroic domains 



provides a versatile and powerful tool for analyzing the 
physical properties of ferroic materials in general, and 
ferroelectrics in particular. Our own analysis of existing 
data suggests that the thickness of 180° domain walls 
in ferroelectrics is extremely narrow (of the order of one 
unit cell), and that regular patterns of sub- micron do- 
mains can be achieved in photonic crystals. More studies 
of domain periodicity as a function of crystal thickness 
should be carried out to establish domain wall thickness 
and exchange parameters for other relevant ferroelectrics 
such as LiNbOs and BaTiO^. 
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